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Job Shop Scheduling 



□ Problem with m machines and n jobs 

□ Each job visits some or all of the machines 

■ Only once (or multiple times if recirculation is allowed) 

□ Each customer order is unique and of small batches 

■ Wafer fabrication in semiconductor industry 

■ Hospital (patients are jobs or machines!) 

□ Special case of project scheduling with workforce 
constraints 

□ Very difficult to solve! (NP-hard) 
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operation/task/ 
activity 



precedence constraint 
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Classic scheduling theory 



□ Look at a specific machine environment with a 
specific objective 

□ Analyze to prove an optimal policy or to show that no 
simple optimal policy exists 

□ Thousands of problems have been studied in detail 
with mathematical proofs! 
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Single machine and parallel machines 



□ Single machine is the simplest models in job shop. 

□ Parallel machine models are equivalent to flexible 
flow shop with a single work center. 

□ These problems can be solved by heuristics 
(dispatching rules). 

□ Sometimes heuristics can be proved to be optimal for 
simple cases 
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Dispatching rules 



□ Dispatch rule can be static or dynamic. 

□ One machine problems (WSPT, EDD, MS, ATC) 

□ Parallel machines (LPT) 

> Prioritize all waiting jobs 

■ job attributes 

■ machine attributes 

■ current time 

> Whenever a machine becomes free: select the job with 
the highest priority 
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Release/due date related 



□ Earliest release date first (ERD) rule 

■ variance in throughput times 

□ Earliest due date first (EDD) rule 

■ maximum lateness 

□ Minimum slack first (MS) rule 

max(<i ; . —p. -t,6) 

■ maximum lateness 

□ Apparent tardiness cost first (ATC) 

■ maximum total weighted lateness 
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□ Longest Processing Time first (LPT) rule 

■ balance load on parallel machines 

■ makespan 

□ Shortest Processing Time first (SPT) rule 

■ sum of completion times 

■ WIP 

□ Weighted Shortest Processing Time first (WSPT) rule 

□ Critical Path (CP) rule 

■ precedence constraints 

■ makespan 
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□ Very simple to implement 

□ Optimal for special cases 

□ Only focus on one objective 

□ Combine several dispatching rules: 

Composite Dispatching Rules 



Joao Miguel da Costa Sousa 



135 



W 



Single machine models & WSPT 



□ n jobs with /? ., r- and dL 

□ Total weighted completion time should be minimized: 

> Solution: Weighted Shortest Processing Time 
(WSPT) first is optimal. 

■ Schedules jobs in decreasing order of w ■//?-. 

□ SPT rule starts with job with the shortest/?-, moves on 
to job with second shortest/?-, and so on. 
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□ Suppose it is not true and schedule S is optimal. 

□ Then there are two adjacent jobs, say job j followed by 
job k such that 



Pj Pk 
□ Do a pairwise interchange to get schedule S' 



t+Pj + Pt 



t+Pj + p k 
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Q Proof 


The weighted 


completion time of the two jobs 
(t+ Pj)Wj + (t + pj + p k )w k 


under S is 


The weighted 


completion time of the two jobs 

(t+P k )w k + {t+p j + p k )w j 


under S ' is 


Then: 






(t + Pj )w j + 


(t + p. + p k )w k =(t + Pj)Wj + pjw k +(t + p k )w k 




>(t + p.)wj + p k Wj + (t + p k )w k 




= (t + p k )w k +(t + p k 


+ Pj)Wj 


Contradicting 


that S is optimal. 
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Single machine models & EDD 



□ = o 



'max 



□ Each job has its own d- 

□ Objective: minimize lateness L r 

□ Earliest Due Date (EDD) results in optimal schedule 
> Order operations in increasing order of df ■ 
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□ WSPT and EDD are static. 

□ Static - basis for ordering operations does not change 
based on scheduling decisions 

■ All operations can be sorted once 

□ Dynamic - scheduling decisions change the order of 
remaining operations 

■ Need to resort operations in queue (potentially) after every 
decision 
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□ = o 

□ Objective: minimize lateness L r 

□ Minimum Slack (MS) orders operations at time t in 
descending order of: 

■ max(df • - Pj - t, 0) 

□ MS does not guarantee optimal schedule! 
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Composite rule 



Two good heuristics: 

□ Weighted Shorted Processing Time (WSPT) 

■ optimal with due dates zero 

□ Minimum Slack (MS) 

■ Optimal when due dates are "spread out" 

□ Any real problem is somewhere in between 

□ Combine the characteristics of these rules into one 
composite dispatching rule. 
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Y Composite Dispatching Rule 



□ One-machine, r = 

□ Objective: minimize weighted tardiness Ew-7^- 

□ Apparent Tardiness Cost (ATC) rule orders 
operations in descending order (K is a parameter): 



/,(*)= 



WSPT 
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Special cases 



□ If K is very large: 

■ ATC reduces to WSPT 

□ If ^T is very small and no overdue jobs: 

■ ATC reduces to MS 

□ If ^T is very small and overdue jobs: 

■ ATC reduces to WSPT applied to overdue jobs 
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Choosing K 



□ Value of K determined empirically 

□ Related to the due date tightness factor 

d 



T = l- 



c 



max 



and the due date range factor 



d —d ■ 

D _ max mm 



c 



max 
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□ Usually 1.5 < K<4.5 

□ Rules of thumb: 

■ Fix K = 2 for single machine or flow shop. 

■ Fix K = 3 for dynamic job shops. 

□ Adjusted to reduce weighted tardiness cost in 
extremely slack or congested job shops 

□ Statistical analysis/empirical experience 
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□ One-machine problem with different r ■ 

□ Objective: minimize lateness L 



'max 



□ Problem is NP-hard 

□ Possible algorithms to solve the problem 

■ Branch-and-bound (see Appendix B of Pinedo's book) 

■ Dynamic programming 
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Parallel machines 



□ A set of m machines in parallel is available. 

□ Objective: minimize makespan C max 

□ Longest Processing Time (LPT) first 

■ pick operations in descending order of processing time 

□ LPT balances the loads of the machines (why?). 

□ LPT does not guarantee optimality. 
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Parallel machines 



> Objective: minimize completion time EC} 

□ SPT assures optimality, even when preemptions are 
allowed. 

> Objective: minimize weighted completion time Ew- Cj 

□ WSPT does not assures optimality. 

> Objective: minimize total weighted tardiness Ew-7^- 

□ This more general problem is even harder. ATC can be 
applied, but solutions can be poor. 
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Dispatching rules for Job Shop 




Q 
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Dispatching rules for Job Shop 




Q 
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, I 


What is the dispatching 
rule being used? 
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Dispatching rules for Job Shop 



JO 



Jl 



M\ — £M 



EQ 



J2 


MJ — a ^E^H 
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5/ Dispatching rules for Job Shop 



JO 
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J2 



M\ 
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Dispatching rules for Job Shop 



bd 



jo El 

" I bJ I — 13 



J2 



[Ml — n 



Q 
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Dispatching rules for Job Shop 



\si 
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J2 



[Ml || 
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□ Job shop with n jobs and m machines. 

□ Each job visits some machines in a given order without 
recirculation. 

□ Processing of job j in machine i is operation (/, j) with 
duration/?^-, and (i,j)eN nodes. 

♦ Objective: minimize makespan C max 

□ Problem can be represented in a disjunctive graph. 
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□ Direct graph G = (N, A, B) with a set of N operations. 

□ Arcs A - conjunctive arcs represent the precedence 
relationships between processing operations of a job. 

□ Arcs B - disjunctive arcs connect two operations 
which belong to two different jobs, that are to be 
processed on the same machine. 

□ Disjunctive arcs form n cliques (in a clique any two 
nodes are connected to one another). 

□ Source U and sink V. 
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Disjunctive graph 



Example of a job shop problem: 4 machines and 3 jobs 
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v Disjunctive graph 



□ Feasible schedule - selection of one disjunctive arc 
from each pair. Each selection of arcs within a clique 
must be acyclic. 

□ Let D be a subset of selected disjunctive arcs. 

□ Makespan of a feasible schedule is the longest path in 
G(D) from the source U to the sink V. 

□ The problem is thus minimizing the longest (critical) 
path. 
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Disjunctive programming 



□ Based on the disjunctive graph. 

□ Let y t j be the starting time of operation (/, j) (operation 
of job j in machine i) 

□ A/"- set of all operations 

□A - set of all conjunctive constraints 

□ 5 - set of all disjunctive constraints 
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v Disjunctive 


programming formulation 


minimize f 

max 








subject to 








^hj yij — Pi] 




for all (i 


J)^(h,j)eA 


C — y.. > p.. 

max y ij r y 




for all (i 


,J)eN 


yu - y ik > Pit or y ik 


— y . > p.- 


for all (i 


,k) and(/, j) 


v.. >0 




for all (i 


J)eN 
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Disjunctive programming formulation 




mi] 
sut 


tiimize r* 

max 

)ject to 




i- ii / • .\ / r -\ 4 


y hj - yu * Py for all (i, j) -> (h, j) e A 




C -y 

max y 

y^ - y ik 

y. >0 


..>/?.. for all (ij)e N 
* Pik or y ik ~ y^ 2. p^ for all (i, k) and (i, j) 

for all (i, j) e N 




An operation cannot start before the 
previous operation (in the job) ends 
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minimize 



subject to 



c 



max 



y^ y^ Pij 



for all (i, j) -> (h, j) e A 



max ^y -^ 



y 



for all (i, j)e N 



yij - y ik > Pik or y tt 

y. >0 

■^y 



y > p 



for all (i,k) and (/, j) 
for all (i, 7) e N 



All operations must end before makespan 
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v Disjunctive programming formulation 




minimize r* 

max 

subject to 
y hj ~ y^ ^ Pij f or all (i, j) -> (h, j) e A 
C^-y,>P, for all (ij)e N 




yy ~ y ik * Pik °r y ik - y i} > P tj for all (/, k) and (/, j) 




y. >0 


for all (i,j)e 


N 




One disjunctive arc must be chosen 
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5/ Disjunctive programming formulation 



C 



minimize 

max 



subject to 

y>v y^ — Py 

C — y.. > p.. 

max J i] f y 

Im - lit > Pik or y ik - yu > Pa 



for all (/,;)-> (/!,;)€ A 

for all (ij) e N 

for all (i,k) and (/, j) 



y >0 



for all (ij)e N 



Start times cannot be negative 
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Solution Methods 



□ Exact solution 

■ Branch & Bound 

■ 20 machines and 20 jobs 

□ Dispatching rules (16+) 

■ Shifting Bottleneck 

□ Search heuristics 

■ Tabu search, Simulated Annealing, Genetic Algorithms, etc. 
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Shifting Bottleneck Heuristic 



□ Minimize makespan in a job shop 

□ Let M denote the set of machines 

□ Let M cMbe machines for which disjunctive arcs 
have been selected 



□ Basic idea 

■ Select a n 

■ Sequence the operations on this machine 



Select a machine in M - M to be included in M 
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Shifting Bottleneck Algorithm 



Step 1: Set the initial conditions 

■ Set M = 0. Graph G is the graph with all the conjunctive 
arcs and no disjunctive arcs. 

■ Set C max (M ) equal to the longest path in graph G. 
Step 2: Analysis of the machines still to be scheduled 

■ For each machine z in M - M : formulate a single machine 
problem with all operations subject to release dates and due 
dates. Release date is the longest path in G from the source 
to the node. Due date is the longest path in G from the node 
to the sink and subtracting p... 

■ Minimize L m „ in each machine. 
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Step 3: Bottleneck selection 

■ The machine with the highest cost is designated the 
bottleneck. 

■ Insert all the corresponding disjunctive arcs in graph G. 

■ Insert machine which is the bottleneck in M . 

Step 4: Re sequencing all machines scheduled earlier 

■ Find the sequence that minimized the cost and insert the 
corresponding disjunctive arcs in graph G. 

Step 5. Stopping condition 

■ If all machines are scheduled (M = M) then STOP, else go 
to Step 2. 
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Example 5.4.2 (p. 89) 



Jobs 


Machines 


Processing times 


1 


1,2,3 


Pll =10,p 21 =8,p 31 =4 


2 


2,1,4,3 


P22= 8 > Pl2= 3 > P42= 5 » P32= 6 


3 


1,2,4 


p 13 =4, p 23 =7, p 43 =3 




Joao Miguel da Costa Sousa 



u 




172 



^ 



Example 5.4.2 (p. 89) 



□ Other form of presenting processing times: 



Ml £] 
M2Q 
M3 
M4 £^ 



Jl 



J2 



J3 
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Initialization 



□ M = set of all machines 

□ M = set of "already scheduled" machines 

■ Initially MO is empty 
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Stepl 



□ Find release date and due date of each operation 

□ Remove all sequence constraints among activities in 
M - M , use CPM to find critical path and minimize 
start time, maximize ending time for each activity 

■ Since M is initially empty, only conjunctive arcs appear in 
the graph. 
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Step 1: Find C max (M ) 














1= 


J — * — L 








l-B 




-3 

C max (M0) 2 
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Jl 




□ 




NOTE: This is a 
job-based Gantt chart 




J2 
J3 








□ 


















Joao 


Miguel da Costa Sousa 








176 







I 


Step 1: Find release and due dates 


o 








C max( M °) 




Jl ► — [ 




[0 10] [10 18] [18 22] 


Jl 


D 




J2 
J3 




j2 — f"H~]~n 


[ m 


[0 8] [8 11] [11 16] [16 22] 

*□- - 

[0 12] [4 19] [11 22] 
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Step 2: Find optimal 1-machine schedules 



□ Using release and due dates, min. L 



max 



L j =C j -d j 
L max = max(L y ) 



Ml I I D D 

[0 10] [8 11] [0 12] 



M2 



Jl 



J2 



J3 



[10 18] [0 8] [4 19] 



M3 



[18 22] [16 22] 



M4 g |j 

[11 16] [1122] 
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11 16 



L ma x(D = 5 



L max (2) = 5 



| L max (3) = 4 

16 22 



L max (4) = 
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Step 3: Add machine to M 



'max 



□ Pick machine with highest L r 

□ Use sequence found in Step 2 

■ Arbitrarily choose to add machine 1 
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Graph after iteration 1 




C max {{ 1 }) = C WM ,(0) + L H ) = 22 + 5 = 27 
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Iteration 2, Step 1 



□ M ={1} 

□ C max (M ) = 27, find release and due dates 



M4 



Ml 




















J2 




J3 




M2 


Jl 




M3 





□ 



Joao Miguel da Costa Sousa 



M2 



Jl 




J2 
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[10 23] [0 10] [17 24] 

M3 

[18 27] [18 27] 

M4 D [ 

[13 21] [24 27] 
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□ Using release and due dates, min. L 



max 



M2 



Jl 



J2 




J3 



[10 23] [0 10] [17 24] 



M3 



[18 27] [18 27] 



M4 



[13 21] [24 27] 



J2 



Jl 


J3 



8 10 18 



L max (2) = 1 
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Iteration 2, Step 3 



□ Pick machine with highest L max 

□ Use sequence found in Step 2 

Anax(2) - Anaxw) = 1 

■ Arbitrarily choose to add machine 2 
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Graph after iteration 2 




10 



\10 



<5> 



/ 3 




C max ({h2}) = C max ({l}) + L max (2) = 27 4- 1 = 28 
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Iteration 3, Step 1 



□ M ={1,2} 

□ C max (M ) = 28, find release and due dates 















Ml 




















M2 


J2 




ji 


J3 




M3 




M4 


n n 
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Iteration 3, Step 2 



□ Using release and due dates, min. L 



max 



M3 



[18 28] [18 28] 



[13 22] [25 27] 
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Iteration 3, Step 2 



□ L max (3) = L max (4) = 

□ Thus, a final schedule was found in Step 2 



Ml 


H 




M2 


J2 


Jl J3 


M3 




M4 


n n 
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□ Very effective 

■ Relatively fast 

■ Good solutions 

■ More general Job Shop problems can be solved as well 

□ 'Just a heuristic ' 

■ No guarantee of optimum 
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Branch-and-Bound 



□ Represent problem as an Integer Programming (IP) 
problem 

■ Sequence of every pair of operations is a 0-1 variable 

□ Use Branch-and-Bound (B&B) to find solution 

□ Will find optimal solution (if given enough time!) 
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v Constraint Programming 

INiTITUTO O CT 



□ B&B (but not IP) plus inference 

□ Every time you branch, use specialized algorithms to 
find other decisions that must be true 

□ May also use sophisticated branching heuristics 

□ May need to backtrack (if the heuristic decision has 
made a mistake) 

□ Also will find optimal if enough time is given. 
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□ On disk with book 

□ Generic job shop scheduling system 

□ User friendly windows environment 

□ C++ object oriented design 

□ Can add own routines 

LI Look at slides from University of Nottingham coming 
with book! or: 

http://www.cs.nott.ac.uk/~sxp/ 
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LEKIN 



P*k;<],1 itMt^i;,^Aiim,,^st 



Workspace File Schedule lob Operation 5 art Tools Window Help 



B*aJ» %\ r ::aBB □ H ffl B 




■ For Heh press Fl 
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